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Abstract: Using the example of the rare decay B — > ^s7, we analyse the importance of 
interference effects for the bounds on the parameters in the squark mass matrices within 
the unconstrained MSSM. In former model-independent analyses no correlations between 
the different sources of flavour violation were taken into account. In our new analysis we 
include the contributions from charged Higgs bosons, charginos and neutralinos and their 
interference effects and, even more important, the effects that appear when several flavour 
violating parameters, i.e. several off-diagonal elements in the squark mass matrices, are 
switched on simultaneously. We derive new bounds on certain off-diagonal elements of 
the squark-mass matrix which are in general one order of magnitude weaker than the 
previous bounds. 
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1. Introduction 

Flavour changing neutral current (FCNC) processes provide crucial guidelines for super- 
symmetry model building. Besides the Cabibbo-Kobayashi-Maskawa (CKM)-induced con- 
tributions, there are generic super symmetric contributions induced by flavour mixing in the 
squark mass matrices within the so-called unconstrained minimal supersymmetric standard 
model (MSSM). The structure of the MSSM does not explain the suppression of FCNC 
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processes, which is observed in experiments; this is the crucial point of the well-known 
supersymmetric flavour problem. 

Among neutral flavour transitions involving the third generation, the rare decay B — > 
X<j7 is at present the most important one |j^], as it is the only inclusive mode that is already 
measured and that already provides theoretically clean and rather stringent constraints 
on the parameter space of various extensions of the SM Although the theoretical SM 
prediction, up to next-to-leading logarithmic (NLL) precision for its branching ratio, is 
in agreement with the experimental data, it is still possible that the rare decay B Xg'j 
leads to the flrst evidence of new physics by a signiflcant deviation from the SM prediction, 
for example in the observables concerning direct CP violation ||5|. 

The decay B — > Xgj is sensitive to the mechanism of supersymmetry breaking because, 
in the limit of exact supersymmetry, the decay rate would be just zero, B{B — > Xg^) = 0. 
Flavour violation thus originates from the interplay between the dynamics of flavour and 
the mechanism of supersymmetry breaking and FCNC processes may contribute to the 
question of which mechanism ultimately breaks the supersymmetry. 

Former analyses in the unconstrained MSSM neglected QCD corrections and only used 
the gluino contribution to saturate the experimental bounds. Moreover, no correlations 
between different sources of flavour violation were taken into account. In this way, one 
arrived at 'order-of-magnitude bounds' on the soft parameters |^ |7[ H- In the sensitivity 
of the bounds on the down squark mass matrix to radiative QCD corrections was analysed, 
including the SM and the gluino contributions. In the new analysis [10| we present here, we 
include the contributions from charged Higgs bosons, charginos and neutralinos and their 
interference effects and, even more important, the effects that result when several flavour 
violating parameters, i.e. several off-diagonal elements in the squark mass matrices, are 
switched on simultaneously. 



2. Phenomenological Analysis 

To understand the sources of flavour violation that may be present in supersymmetric mod- 
els, in addition to those enclosed in the CKM matrix, one has to consider the contributions 
to the squark mass matrices 

2 _ ( /, LL + FfLL + Df LL (m^ ^ lr) + LR 

\ y^f^LR) +F}RL m^f,RR + FfRR + DfjiR 

where / stands for up- or down-type squarks. We recall that the matrices tUu^ll and 
f^d^LL cannot be specified independently; SU{2)l gauge invariance implies that mu,LL = 
Km^ LLK^ where K is the CKM matrix. In the super-CKM basis, where the quark mass 
matrices are diagonal and the squarks are rotated in parallel to their superpartners, the F 
terms from the superpotential and the D terms turn out to be diagonal 3x3 submatrices 
of the 6x6 mass matrices This is in general not true for the additional terms rrij, 

originating from the soft supersymmetry breaking potential. Because all neutral gaugino 
couplings are flavour diagonal in the super CKM basis, the gluino contributions to the decay 
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7 are induced by the off-diagonal elements of the soft terms fn^f i^j^i ^"f^jiji^ ^'/^rl' 
Since there are different contributions to this decay, with different numerical impact on 
its rate, some of these flavour violating terms may turn out to be poorly constrained. 
Thus, given the generality of such a calculation, it is convenient to rely on the mass 
eigenstate formalism, which remains valid even when some of the intergenerational mixing 
elements are large, and not to use the approximate mass insertion method (MIA), where 
the off-diagonal squark mass matrix elements are taken to be small and their higher powers 
neglected. In the latter approach the reliability of the approximation can be checked only 
a posteriori. 

As a first step, it is convenient to select one possible source of flavour violation in the 
squark sector at a time and assume that all the remaining ones are vanishing. It should be 
stressed that one already excludes any kind of interference effects between different sources 
of flavour violation in this way. Following ref. Q, all diagonal entries in j^j^, m?^ 
and m\ are set equal and their common value is denoted by m?. The branching ratio 
can then be studied function of 

m| m| m| 

We find that only those parameters get stringently bounded by i? — > ^^7, which can 
generate contributions to the two five-dimensional gluino-induced dipole operators, namely 
Sd,LR,23 and 6d,RL,23- The two corresponding operators are connected by chirality and are 
denoted by 07^,3 and O-rg^g in the following. As the gluino yields, intrinsically, the dominant 
contribution by far, we also find that the bounds on 6(1,lr,23 and 5d,RL,23 are only marginally 
modified by chargino, neutralino and charged Higgs boson contributions. 

In the second part of our analysis, we investigate whether the obtained bounds re- 
main stable if all off-diagonal elements, which induce the decay B — > ^^7, are varied 
simultaneously. For our analysis we explore various scenarios that are characterized by 
the values of the parameters Mj;^- , tan/3, Mgusy, ''T^g. We regard this as reasonable, be- 
cause we expect that these input parameters, which are unrelated to fiavour physics, will 
be fixed from flavour conserving observables in the next generations of high energy ex- 
periments (provided low energy SUSY exists). We note that the common SUSY scale, 
Mgusy, flxes in our scenarios the general soft squark mass scale rriq (see (p.2[) ) and the 
first diagonal element of the chargino mass matrix M2. The mass of the charged Higgs 
boson Mij~ are fixed to be Mjj- = 300 GeU. We also allow for a non-degeneracy of 
the diagonal elements in the matrices m\ j^j^, and j^j^. To implement this, 

we define (5-quantities in addition to those in eq. ( |2.2[ ), which parametrize this non- 
degeneracy: 5f^LL,n = {{rn^f^LLlii ~ ^ud analogously for df^RR^u. In our Monte 
Carlo analysis these diagonal (5-parameters are varied in the interval [—0.2,0.2]. On the 
other hand, the off-diagonal ones (in eq. ( |2.2| )) are varied in the interval [—0.5,0.5], by 
use of a Monte Carlo program. There are, however, two exceptions. First, we do not 
vary those off-diagonal 5's with an index 1; the latter (5's we set to zero, since they are 
severely constrained by kaon decays (see for example [^). Second, also ("i^^^)33 is not 
varied, but fixed such that the mass of the lightest neutral Higgs boson gets heavy enough 
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Figure 1: Contours in the Sd.LR,23-Sd,RL,23 plane. In tlie left frame, Sd.LR.23 and Sd,RL,23 are the 
only flavour violating parameters. In the right frame, Sd,LR.23, 5d.RL.23, Sd.LL.,23 and Sd.RR.23 are 
flavour violating parameters; Sd.LR,22, Sd,LR,33 are also non-zero. We only consider SM and gluino 
contributions and the other parameters are fi = 300 GeV, tan/3 — 10, Afgusy — 500 GeV, x = 1. 



to be compatible with experimental bounds [11|. We plot those events, corresponding to 
2.0 X 10""^ < BR{B Xsj) < 4.5 x 10"^, which is the range allowed by the CLEO 
measurement. We start with the following parameter set: /i = 300 GeV, Mjj- = 300 GeV, 
tan/3 = 10, Msusy = 500 GeV, x = / M^^^^ = 1 and Xt = 750 GeV. In fig. |, we 
only consider SM and gluino contributions. In the left frame we present the constraints on 
Sd,LR,23 and 5d,RL,23 when these are the only flavour-violating soft parameters; the diagonal 
(5-parameters defined above are also switched off. As expected, stringent constraints are 
obtained. The hole inside the dotted area represents values of 5rf,L_R,23 and 5d^RL,23 for 
which the branching ratio is too small to be compatible with the measurements. In the 
right frame we now investigate interference effects from different sources of flavour viola- 
tion, where we allow for non-zero 5d,LR,23, Sd,RL,23, Sd,LL,23, Sd,RR,23, Sd,LR,22, and 5d,LR,33- 
All these parameters are varied between ±0.5. As can be seen, the bounds on 5d,LR,23 and 
Sd,RL,23 get destroyed dramatically (note the different scale in the two plots). The reason 
for this is that there are now new contributions to the five-dimensional dipole operators. 
As an example, the combined effect of 5d,LR,33 and 5d,LL,23 leads to a contribution to the 
Wilson coefficient of the gluino- induced magnetic dipole operator Ojg^g. The sign of this 
contribution can be different from the one generated by 5d^LR,23- As a consequence, the 
bound on 5d^LR,23 gets weakened. To illustrate this more quantitatively, we assume for the 
moment that there are only these two sources that can generate Ojg^g, i.e. we switch off 
the other (5-quantities. If 5d,LR,23 is larger than the individual bound from the first part of 
the analysis, it is necessary that the product of 5d,LR,33 and 5d,LL,23 also be relatively large; 
only in this case can the two sources lead to a branching ratio compatible with experiment. 
This feature is illustrated in fig. ^; only values of 5d,LR,23 and values of 5d,LR,33 • Sd,LL,23 
that are strongly correlated lead to an acceptable branching ratio. As clearly visible from 
fig. |2|, the correlation between the two sources for Ojg^g is essentially linear. This implies 
that the linear combination com := 5d,LR,23 + / ^d,LR,33 ■ ^d,LL,23 gets constrained if / 
is chosen appropriately. Stated differently, the Wilson coefficient of the five-dimensional 
magnetic dipole operator is essentially proportional to the combination com defined above. 
This implies in turn, that for the values of the parameters we are using at the moment, 
the Wilson coefficient is well approximated by its double mass insertion expression. Thus, 
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Figure 2: The parameters Sd.LR.23 and Sd.LR.ss • Sd.LL.23, which are compatible with the data on 
B Xsj, are shown by dots. Values lying on the solid line lead to a vanishing contribution of the 
five-dimensional magnetic dipole operator Or^.g in the MIA. See text. 

the coefficient / can be fixed analytically (see |jl^); the numerical values of f{x) for some 
values of x read 0.74 for x = 0.3, 0.68 for x = 0.5, 0.60 for x = 1.0 and 0.52 for x = 2.0, 
respectively. The solid line in fig. |2| represents pairs {6d,LR,23, ^d,LR,33 ' ^d,LL,23) for which 
the combination com is zero. The points scattered around this line therefore represent 
Monte Carlo events for which this combination is small. We now turn back to the scenario 
of fig. |, in which all the parameters 5d,LR,23, h,RL,23i Sd,LL,23, ^d,RR,23, ^d,LR,22, ^d,LR,33 
are varied simultaneously. In this case, the linear combinations 

LCi = 5ci,RL,23 + f{x)^d,RR,23 ' ^d,RL,33 + f{x)^d,RL,22 ' ^d,LL,23i 

LC2 = Sd,LR,23 + f(.x)5d,LR,22 ' ^d,RR,23 + f{x)5d,LL,23 ' ^d,LR,33^ (2-3) 

are expected to get constrained. In fig. ^ we show the allowed region for LCi and LC2. 
There, we allow all non-diagonal 5-parameters to vary between ±0.5. In addition, we also 
allow for non-equal diagonal soft entries, by varying the parameters and Sf^RR^u 

between ±0.2. With the latter choice we still guarantee the hierarchy between diagonal 
and off-diagonal entries, but we get rid of the unnatural assumption of degenerate diagonal 
entries. In the left frame, we include only the SM and gluino contributions. We find that 
the linear combinations LCi and LC2 indeed get stringently bounded. In the right frame of 
fig. ^ we test the resistance of these bounds when the additional contributions (i.e., those 
from charginos, charged Higgs bosons and neutralinos) are turned on. In this case also 
Su,LR,23, Su,RL,23, 6u,LL,23, Su,RR,23 and 5u,LR,22 are varied in the range ±0.5. We find that 
the bound on LCi remains unchanged, while the one on LC2 gets somewhat weakened. 
This feature is expected, because charginos and charged Higgs bosons contribute to the 
unprimed operator. 

At this point we should stress that these plots were obtained by choosing the renormal- 
ization scale = 4.8 GeV and by requiring all squark masses to be larger than 150 GeV. 
We checked that the bounds on LCi and LC2 remain practically unchanged when the 
renormalization scale is varied between 2AGeV and 9.6 Gel^; they are also insensitive to 
the value of the required minimal squark mass, as we found by changing msquarkmin from 
150 Gey to 100 Gey or 250 GeT^. Moreover, we also checked whether the restriction to the 
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Figure 3: Contours in tlie LC1-LC2 plane with (5d,z,ij,23, Sd,RL,23, ^d,LL,23, ^d,RR,2i, Sd,LR,22, 

Sd,LR,3a, SuM.L,2a, Su,LL,23, Km.r,23, and Su^LR,22 all non- vanishing. In the left frame, we 

consider only SM and gluino contributions whereas in the right frame we also include chargino, 
charged Higgs boson and neutralino contributions. 

H = +300 GeV scenario is too severe: we redid the complete analysis for n = —300 GeV 
and confirmed that there are no differences between the results with these two choices. 

There is a remark in order. If we got rid of the hierarchy of diagonal and off-diagonal 
entries in the squark mass matrices, stringent bounds on the simple combinations LCi and 
LC2 certainly would no longer exist, simply because there would then be more contributions 
to the five-dimensional operators of similar magnitude. In this case, however, the full 
Wilson coefficients of the five-dimensional operators would still be stringently constrained 
by the experimental data on S — > Xsj. Unfortunately, in this case not much information 
can be extracted for the individual soft parameters or simple combinations thereof. 

Finally, we extend our analysis to other values of the input parameters. Wc analyse 
the bounds on the soft parameters within the following parameter sets: Mgusy = 300 GeV^, 
500 GeV, 1000 GeF. For tan/3 we explore the values: tan/3 = 2, 10, 30, 50. Furthermore, 
the gluino mass rUg is varied over the values x = rn^ / M'^^^y = 0.3 , 0.5 ,1,2. Surprisingly, 
we find that the constraints on LC\ and LC2 are completely stable over large parts of 
the parameter space. However, the bounds get weakened when tan /3 values as large as 50 
are chosen. This effect gets enhanced when the general mass scale rUq in the squark mass 
matrices decreases with the parameter Mgusy There are two main reasons why the bounds 
get weakened in these scenarios. First, in the large tan/3 regime, the term (-Fd,Li?)33 gets 
strongly enhanced because it is proportional to tan/3. Particularly, for tan/3 = 50 and 
-^susy = 300 GeF, the term (F(i,Li?)33 is of the same magnitude as the diagonal entries 
of the squark mass matrix. Thus, the contributions to the Wilson coefficients of the five- 
dimensional gluino operators (induced by (i^d,LR)33 in combination with 5d,LL,23 or S(1,rr,23) 
become important enough to weaken the bounds on LGi and LC2 significantly. The relative 
importance of this F term is of course increased if the general soft squark mass scale Msusy 
is decreased. Second, within the large tan/3 regime the contributions from charginos get 
enhanced and therefore also weaken the bounds on LG2. Because the parameter (Fd,Lii)33 
is actually proportional to the product of tan/3 and we conclude from the above findings 
that the bound on LCi is unchanged if we increase the value of fi and decrease the value of 
tan (3 so that the product of both parameters is constant; the bound on LC2 is then even 
stronger, because the chargino contribution is smaller for increasing /i. 
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3. Conclusions 



Our new model-independent analysis of the rare decay B Xsj, based on a system- 
atic leading logarithmic QCD analysis, mainly explored the interplay between the various 
sources of flavour violation and the interference effects of SM, gluino, chargino, neutralino, 
and charged Higgs boson contributions. In former analyses, no correlations between the 
different sources of flavour violation were taken into account. Unlike previous work, which 
used the mass insertion approximation, we used in our analysis the mass eigenstate for- 
malism, which remains valid even when some of the intergenerational mixing elements are 
large. We singled out two simple combinations of elements of the soft parts of the down 
squark mass matrices, which stay stringently bounded over large parts of the supersymmet- 
ric parameter space, excluding the large tan /3 and the large /i regime. These new bounds 
are in general one order of magnitude weaker than the bound on the single off-diagonal 
element 6d^LR,23j which was derived in previous work 12| by neglecting any kind of in- 
terference effects. It seems that the flavour problem is less severe in the B system than 
often stated. Thus, it would be interesting to explore also the consequences of natural 
interference effects between different sources of flavour changes within the kaon sector. 
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